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Analytical  results  of  rotor  blade-vortex  interaction  airloads  are  presented  with  two  different  wake  models  in 
the  eomprehensive  analysis  CAMRAD  II,  and  these  ealculated  results  are  eompared  with  the  experimental  data 
obtained  from  the  higher-harmonic-control  aeroacoustic  rotor  test  (HART-I)  program.  The  HART  rotor  was 
a  40%,  Mach-scaled  model  of  the  hingeless  BO-105  main  rotor.  Two  wake  models  used  in  the  eomprehensive 
analysis  are  the  single-trailer  and  the  multiple-trailer  free-wake  models.  The  multiple-trailer  wake  model  shows 
good  predietion  of  lift  distribution,  M^c„,  as  a  function  of  azimuth  for  baseline,  minimum  noise,  and  minimum 
vibration  eases,  and  shows  signilieant  improvement  relative  to  the  single-trailer  model  in  predietion  of  the  tip 
vortex  wake  geometry. 


Nomenclature 


= 

blade  section  normal  force  coefficient 

deoils 

= 

fraction  rollup 

G 

total  strength  of  trailed  vorticity 

''^cons 

consolidation  scaling  factor 

M 

= 

Mach  number 

r 

= 

blade  radial  station 

rc 

centroid  of  vorticity 

ro 

moment  (radius  of  gyration)  of  vorticity 

ois 

rotor  shaft  angle  (positive  aft) 

r 

= 

blade  bound  circulation 

T 

= 

wake  age  for  the  trailed  vortex  filament 

^cons 

consolidation  time  constant 

Introduction 

LADE-VORTEX  interactions  (BVI)  are  the  source  of  an¬ 
noying  and  intrusive  noise  from  helicopter  rotors.  This  phe¬ 
nomenon  is  one  of  the  distinctive  features  of  rotor  wakes,  and  the 
noise  source  becomes  dominant  during  low- speed  descent  and  ma¬ 
neuvering  flights,  where  the  rotor  wake  is  blown  back  into  the  rotor 
plane.  High  noise  levels  produced  by  BVI  may  prevent  rotorcraft 
from  achieving  wide  acceptance  for  civil  applications. 

In  the  last  decade,  many  investigations  have  been  conducted  to 
demonstrate  the  potentially  high  payoff  in  noise  and  vibration  re¬ 
duction  using  higher-harmonic  blade  pitch  control  (HHC)  inputs. 
Although  the  HHC  technique  was  originally  intended  for  vibration 
reduction,  several  wind-tunneland  flight  tests  have  shown  BVI  noise 
reduction  of  up  to  10  dB  with  the  application  of  HHC  inputs,  while 
vibration  and  low-frequency  noise  levels  were  increased.'"^  How¬ 
ever,  with  certain  HHC  schedules  (amplitude  and  phase),  simultane¬ 
ous  reductionof  noise  and  vibrationhas  been  found  to  be  achievable. 
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Better  physical  understanding  and  analytical  models  are  necessary 
to  explain  the  basic  physics  of  the  noise  and  vibration  reduction. 

In  1994,  researchers  from  DLR,  German  Aerospace  Research 
Center,  French  ONERA,  NASA  Langley  Research  Center, 
the  German-Netherlands  Wind  Tunnel  (DNW),  and  the  U.S. 
Army  Aeroflightdynamics  Directorate  conducted  the  first  higher- 
harmonic  aeroacoustics  rotor  test  (HART)  rotor  test  with  a  40% 
Mach-scaled  model  of  a  BO- 105  main  rotor  in  the  open-jet  ane- 
choic  test  section  of  the  DNW.^“"  In  this  HART-I  test,  extensive 
measurements  of  acoustics,blade  pressures,  and  blade  deformations 
were  made,  as  well  as  limited  rotor  wake  measurements  each  at  one 
azimuth  angle  on  the  advancing  side  and  the  retreating  side.  More 
detailed  information  about  the  measurements  is  available  in  Ref.  8. 

Previous  correlation  efforts  with  the  HART-I  data,'^“'^  revealed 
that  accurate  prediction  of  2-per-revolution  harmonic  component 
of  the  lift  distribution  is  a  very  challenging  task.  Some  adjust¬ 
ments  to  comprehensive  analysis  models  were  necessary  to  achieve 
reasonable  accuracy  of  the  prediction,  such  as  prescribing  blade 
torsion  from  the  measured  data  or  vortex  wake  geometry  iterative 
refinement. 

The  purpose  of  this  paper  is  to  explore  a  multiple-trailer  free- 
wake  model  in  CAMRAD  II*^  '^  for  further  improvement  of  the 
rotor  airload  prediction  capability. 

Experimental  Data 

Three  HART-I  cases  selected  for  comparison  with  the  experimen¬ 
tal  data  were  the  baseline  case  [(BL)  run  1 40] ,  minimum  noise  case 
[(MN)run  1 38],  and  minimum  vibrationcase  [(MV)  run  133]. These 
cases  were  for  the  same  descent  flight  condition  with  an  advance 
ratio  of  0.15,  a  shaft  angle  of  5.3-degaft  (4.2-deg  aft  when  adjusted 
for  wind-tunnel  correction),  a  hover  tip  Mach  number  of  0.64,  and  a 
thrust  coefficient  of  0.0044.  The  BL  case  is  a  normal  descent  flight 
condition  without  higher-harmonic  pitch  control  inputs.  The  MN 
case  has  3-per-revolutionHHC  inputs  with  0.87-deg  amplitude  and 
296-deg  phase  shift,  which  generates  lower  noise.  The  MV  case  has 
3-per-revolutionHHC  inputs  of  0.83-deg  amplitude  and  178-deg 
phase  shift,  which  produces  MV.  More  detailed  information  about 
the  test  data  is  available  in  Ref.  8. 

CAMRAD  II  Rotor  Modeling 

The  rotor  blade  experiences  strong  induced  velocity  fields  due  to 
the  close  passage  of  tip  vortices,  which  have  been  trailed  from  the 
blade  itself  or  the  other  blades.  For  very  close  interactions,  viscous 
effects  can  be  important,  and  transonic  effects  can  be  also  important 
at  high-speed  flights.  These  effects  could  be  modeled  by  the  use  of 
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sophisticated  methods  such  as  fuii-biown  computationai  fluid  dy¬ 
namics  (CFD)  methods,  but  due  to  computationai  compiexity  and 
high  computation  time  requirement,  the  usage  of  these  methods  has 
been  iimited  for  rotorcraft  airioads  caicuiation.  Instead,  the  iift-iine 
theory  has  been  frequentiy  used  for  generai  rotor  appiication,  ac¬ 
companied  with  a  compiex  rotor  wake  modei  and  with  various  cor¬ 
rections  to  reflect  compiex  physics  of  unsteady  rotor  aerodynamics, 
which  has  been  typicaiiy  a  part  of  the  comprehensive  rotorcraft 
anaiysis. 

For  a  rotor  structure,  the  biade  was  structuraiiydiscretizedinto  16 
nonlinearbeam  elements  (flnite  element),  with  more  beam  elements 
aligned  inboard  to  capture  the  blade  structural  behavior  accurately. 
The  beam  element  was  representedby  3  translational(axial,  lead  lag, 
and  flap)  and  3  correspondingrotational  degrees  of  freedom  (DOF) 
andwouldgenerallyresultin  15  DOFs  foreachbeam element.  For  an 
aerodynamics  model,  a  reflned  aerodynamic  model  was  employed 
using  16  aerodynamic  panels  on  each  blade  to  capture  BVI  more 
accurately,  and  aerodynamic  panels  were  more  densely  distributed 
outboard. 

Unsteady  lift  and  moment  in  the  attached  flow  are  calculated 
based  on  compressible  thin-airfoil  theory.  Because  the  steady  air¬ 
loads  are  obtained  from  airfoil  tables  and  the  lift  deflciency  effects 
are  usually  accounted  for  in  the  wake-induced  velocity,  only  the 
unsteady,  noncirculatory  effects  need  to  be  included  for  a  total  air¬ 
load  calculation.  The  unsteady  aerodynamics  model  selected  in  the 
CAMRAD  II  code  was  based  on  the  modifled  ONERA  EDLIN 
theory.'^  The  modiflcations  were  mainly  for  less  complexity  when 
interfaced  with  the  comprehensive  analysis  code,  such  as  an  ex¬ 
clusion  of  steady  loads  calculation,  a  calculation  about  the  quarter 
chord  for  unsteady  airloads,  or  an  extension  to  the  reverse  flow  re¬ 
gion.  The  lift  effect  produced  by  the  steady,  quarter  chord  upwash 
was  also  added  as  a  part  of  the  unsteady  airloads  calculation. 

The  rotor  vortex  filaments  are  convected  from  each  rotating 
blades,  and  the  wake  is  dominated  by  the  blade  tip  vortices  that 
are  trailed  from  the  tips  of  each  rotor  blade.  The  corresponding  ro¬ 
tor  wake  trajectories  are  described  in  terms  of  the  blade  rotational 
speed  as  well  as  the  free  stream  velocity.  The  induced  velocity  is  cal¬ 
culated  from  the  Biot-Savart  law,  by  integrating  over  the  complete 
wake  surface, 

r  f  X  da 

An  |rP 

where  the  vector  f  is  at  a  distance,  dcr,  from  the  vortex  filament  to 
the  aerodynamic  computation  point  and  T  is  the  vortex  strength  of 
the  trailed  filaments. 

For  the  rotor  wake  structure,  the  trailed  vortices  are  modeled  by 
two  different  free-wake  approaches  in  CAMRAD  II:  a  single-trailer 
wake  model  and  a  multiple-trailer  wake  model.  The  single-trailer 
free-wake  model  includes  a  rolled-up  vortex  model  that  is  assumed 
to  be  formed  at  the  tip  of  the  blade.'*  The  multiple-trailerfree-wake 
model  has  discrete  trailed  vortex  filaments  convected  from  each  of 
the  aerodynamic  panel  edges.  These  trailed  filaments  at  the  panel 
edges  can  be  consolidated  into  a  single  rolled-up  line,  using  the 
trailed  vorticity  moment  to  scale  the  rate  of  rollup. 

The  basic  assumption  of  the  consolidationmodel  is  that  adjacent 
trailed  vortex  lines  of  the  same  sign  eventually  roll  up  into  a  single 
vortex,  located  at  the  centroid  of  the  original  vorticity  distribution. 
The  trailed  vorticity  is  partitioned  into  sets  of  adjacent  lines  that 
have  the  same  sign.  The  total  strength  G,  centroid  and  moment 
(radius  of  gyration)  rg  of  the  trailed  vorticity  in  the  set  are 


,  f  ar 

G4  =  /  -—■(r-rc)^dr 


where  T  is  the  blade  bound  circulation.  The  characteristictime  r^/G 
is  taken  as  a  measure  of  the  rate  of  consolidation.  Thus,  for  trailed 
vortex  elements  at  wake  age  t,  the  total  strength  and  moment  are 


evaluated  at  the  time  (t  —  r)  that  the  vorticity  is  created  and  are  used 
to  calculate  the  time  constant 

Aons  ^cons  ‘  j  G 

Then  the  fraction  of  rollup  is 

/cons  (U  “  Tg)  j 

The  consolidation  begins  at  wake  age  x^.  The  consolidationcan  be 
accomplished  by  entrainment  or  by  compression.  With  the  entrain¬ 
ment  form,  vortex  lines  are  consolidatedinto  a  single  line  of  strength 
/consG-  With  the  compression  form,  /cons  is  the  fraction  of  consol¬ 
idation.  For  the  results  presented  here,  the  entrainment  form  was 
used,  with  fccons  =  1.0  and  n  =  0.33.  Options  to  evaluate  /cons  using 
exponential'^-^"  and  linear  (n  =1)  functions  were  also  available, but 
not  used  here. 

Results  and  Discussion 

The  HART-I  rotor  was  modeled  using  the  CAMRAD  II  compre¬ 
hensive  code.'"-'’  The  rotor  is  a  40%,  Mach-scaled  model  of  the 
hingeless  BO- 105  main  rotor  with  a  radius  of  2  m  and  a  root  cutout 
of  0.35  m,  operating  at  a  nominal  speed  of  1040  rpm.  The  rotor 
blade  has  a  standard  rectangular  tip  with  a  solidity  of  0.077.  It  has 
a  modifled  NACA  23012  airfoil  with  a  chord  length  of  0. 121  m  and 
—8  deg  of  linear  twist. 

Calculations  for  three  test  cases,  BE  (run  140),  MN  (run  138), 
and  MV  (run  133),  were  made.  Eigure  1  shows  the  wake  geome¬ 
try  obtained  using  the  multiple-trailer  model  (entrainment  form). 
Figure  1  shows  the  wake  geometry  behind  each  of  the  four  rotor 
blades,  with  the  reference  blade  at  1 10-deg  azimuth.  The  consolida¬ 
tion  distance  of  the  vortex  filaments  is  longer  behind  the  advancing 
side  than  behind  the  retreating  blade,  which  reflects  a  large  value  of 
the  consolidation  characteristic  time  /  T  on  the  advancing  side. 
Unlike  the  single-trailer  wake  model,  it  is  found  that  vortex  filament 
bundling  in  the  multiple-trailer  vortex  model  occurs  not  only  near 
the  blade  tip  but  also  in  the  blade  midspan.  Recall  that  adjacent 
trailed  vortex  lines  with  the  same  sign  are  consolidatedinto  a  single 
rolled-up  trailer.  Thus,  when  the  loading  on  the  blade  is  negative, 
there  are  three  rolled-up  trailers  in  the  consolidated  wake. 

For  the  trimmed  solution,  the  rotor  was  trimmed  to  the  thrust  and 
the  hub  aerodynamic  moments  (pitching  and  rolling).  The  magni¬ 
tudes  of  the  measured  hub  moments  were  small.  For  the  MN  and 
MV  cases,  the  3 -per-re volution  pitch  control  input  was  added,  and 
then  the  rotor  was  retrimmed  to  reach  the  trim  targets. 


-2-10  1  2 


X(m) 

Fig.  1  Multiple-trailer  wake  geometry  eonvected  from  eaeh  of  the 
rotor  blades  in  the  entrainment  form  (reference  blade  at  20  deg  of 
azimuth). 
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Blade  Section  Airloads 

Blade  section  airloads  were  obtained  by  integrating  the  measured 
blade  surface  pressures  along  the  airfoil  chord  at  the  fully  instru¬ 
mented  radial  stations.  For  the  BL  case,  the  measured  nondimen- 
sional  normal  forces  M^c„  are  shown  as  functions  of  azimuth  angle 
in  Fig.  2,  at  the  blade  span  locations  of  87  and  97%.  The  measured 
data  show  strong  2-per-revolution  variation  of  lift  on  the  front  of  the 
rotor  disk  and  sharp  fluctuations  of  the  BVI  induced  loading  on  the 
advancing  and  retreating  sides.  Calculated  lift  distributions  using 
the  single-trailer  and  multiple-trailer  free-wake  models  are  plot¬ 
ted  with  the  measured  data.  The  CAMRAD  II  single-trailer  wake 
model  producespoorresults  for  the  2-per-revolutionlift  distribution. 
Although  the  CAMRAD  II  calculations  are  slightly  better  outboard 
(97  %),  there  is  a  significant  discrepancy  in  the  2-per-revolutionload- 
ing.  A  similar  trend  was  observed  in  earlier  2GCHAS  calculation 
efforts.'^  ''*  The  multiple-trailer  free-wake  model,  however,  shows 
good  lift  prediction  for  both  the  span  locations.  These  results  imply 


that  the  discrepancy  observed  in  HART-I  airload  correlation  efforts 
is  caused  by  the  use  of  single-trailer  wake  models  in  the  compre¬ 
hensive  codes. 

For  the  MN  case,  the  calculated  results  are  compared  with  the 
measured  nondimensional normal  force  in  Fig.  3,  at  the  blade  span 
locationsof  87  and  97%.  The  measured  data  show  a  much  larger  third 
harmonic  variation  of  the  lift  compared  with  the  baseline  case,  and  in 
fact,  the  peak-to-peakmagnitude  in  the  minimum  noise  case  is  more 
than  doubled  from  the  baseline  magnitude.  This  large  increase  in  the 
peak-to-peakmagnitude  presumably  is  mainly  caused  by  the  large 
3-per-revolutionpitch  control  input.  The  calculated  lift  shows  good 
waveform  prediction,  but  neither  wake  model  was  able  to  match  well 
the  peak-to-peak  amplitude.  The  multiple-trailer  free-wake  model 
gives  a  slight  improvement  of  lift  correlation,  compared  with  the 
single-trailerwake  model. 

For  the  MV  case,  shown  in  Fig.  4,  the  trends  are  similar  to  those 
found  in  the  MN  case.  Much  larger  higher  harmonics  of  lift  are 


Azimuth,  deg 


Azimuth,  deg 


a)  HART  I,  BL,  rlR  =  0.87 


b)  HART  I,  BL,  rlR  =  0.97 


Fig.  2  Calculation  of  M^c„  for  BL  case  with  single-trailer  and  multiple-trailer  wake  models  at  advance  ratio  0.15. 


a)  HART  I,  MN,  r/R  =  0.87  b)  HART  I,  MN,  r/R  =  0.97 

Fig.  3  Calculation  of  M^c„  for  MN  case  with  single-trailer  and  multiple-trailer  wake  models  at  advance  ratio  0.15. 


Azimuth,  deg 


a)  HART  1,  MV,  r/R  =  0.87  b)  HART  1,  MV,  r/R  =  0.97 

Fig.  4  Calculation  of  M^c„  for  MV  case  with  single-trailer  and  multiple-trailer  wake  models  at  advance  ratio  0.15. 


1126 


LIM,  YU,  AND  JOHNSON 


seen  in  the  measured  data,  and  the  calculated  peak- to-peak  magni¬ 
tude  is  well  below  the  measured  data.  This  source  of  the  difference 
between  the  measured  and  calculated  airloads  for  the  MN  and  MV 
is  not  known. 

Blade  Deflection 

Figure  5a  shows  correlation  of  the  lead-lag  deflection  at  the  blade 
tip  for  the  three  test  cases.  The  CAMRAD  II  calculations  are  made 


utilizing  single-trailer  and  multiple-trailer  free-wake  models.  The 
waveform  of  the  lead-lag  deflection  is  well  predicted  for  all  three 
cases,  but  the  calculated  mean  deflection  is  shifted  toward  the  lead¬ 
ing  edge  by  more  than  one  quarter-chord  length  (1  quarter  chord  = 
3  cm).  The  two  wake  models  produce  essentially  the  same  results  for 
the  lag  deflection.  The  difference  between  measured  and  calculated 
mean  lag  deflection  is  probably  produced  by  an  uncertainty  in  the 
blade  torque  offset  or  chordwise  center-of-gravity  position. 


HART  I,  Baseline 


HART  I,  Minimum  Noise 


HART  I,  Minimum  Vibration 


Azimuth,  deg 

Azimuth,  deg 

Azimuth,  deg 

a)  Lead-lag  deflection  at  blade  tip 

HART  1,  Baseline 

HART  1.  Minimum  Noise 

HART  1,  Minimum  Vibration 
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b)  Flap  deflection  at  blade  tip 
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d)  Torsion  deflection  at  60-deg  azimuth  (as  function  of  blade  span) 

Fig.  5  Calculated  blade  response  for  BL,  MN,  and  MV  cases,  using  single-trailer  and  multiple-trailer  wake  models,  at  an  advance  ratio  of  0.15. 
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The  flap  deflection  at  the  blade  tip  is  shown  in  Fig.  5b  for  the  three 
cases  (BL,  MN,  and  MV).  The  peak-to-peak  amplitude  and  phase 
of  the  flap  deflection  are  well  predicted.  As  observed  similarly  in 
prediction  for  the  lag  deflection,  the  two  free- wake  models  produce 
essentially  the  same  results  for  the  flap  deflection. 

In  earlier  studies,'''-'^  the  elastic  torsion  was  found  to  be  critical 
to  accurate  airload  prediction,  whereas  lead-lag  and  flap  deflections 
were  less  important.  Figures  5c  and  5d  show  the  elastic  torsion  at  the 
blade  tip  as  a  function  of  azimuth  angle  and  blade  span,  respectively. 
The  calculated  mean  values  of  elastic  torsion  are  low  for  all  three 
(BL,  MN,  and  MV)  cases,  and  also  the  peak-to-peak  magnitudes, 
are  smaller  than  measured,  for  both  the  multiple-trailer  and  the 
single-trailer  wake  models.  The  spanwise  variation  of  the  elastic 
torsion  is  almost  linear,  so  that  the  difference  between  the  measured 
and  calculated  elastic  torsion  at  the  blade  tip  is  fractionally  spread 
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along  the  span.  The  azimuthal  differences  between  measured  and 
calculatedelastictorsion  are  consistentwith  the  differencesbetween 
the  section  airloads,  shown  in  Fig.  2.  The  measured  torsion  deflec¬ 
tion  exhibits  more  2-per-revolution  content  for  the  BL  case,  and 
more  3 -per- re  volution  content  for  the  MN  and  MV  cases.  Because 
the  3-per-revolution  response  in  the  MN  and  MV  cases  is  caused 
by  the  3 -per- re  volution  pitch  input,  there  probably  is  some  feature 
of  the  blade  structural  dynamics  that  has  not  been  properly  modeled 
in  the  calculations,  which  leads  to  the  underpredictionof  the  torsion 
and  airloads. 

Geometry  of  Wake  Tip  Vortex 

Measurements  of  selectedtip  vortexwake  segments  were  reported 
at  the  azimuth  angles  of  35  and  295  deg  in  Ref.  10.  For  the  BL  case. 
Fig.  6  shows  correlation  of  tip  vortex  wake  geometry  in  the  top 
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b)  Retreating  side,  ^  =  295  deg 

Fig.  6  Comparison  of  measured  and  calculated  tip  vortex  filament  locations  with  single-trailer  and  multiple-trailer  wake  models  for  BL  case. 
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view  and  side  view,  using  both  single-trailer  wake  and  multiple- 
trailer  wake  models  in  CAMRAD  II.  The  top  view  indicates  the 
location  of  tip  vortices  where  significant  BVIs  are  considered  to 
occur,  and  the  side  view  determines  the  miss  distance  of  the  tip  vor¬ 
tices.  The  miss  distance  is  the  vertical  distance  between  the  passing 
blade  and  tip  vortices.  The  fifth  and  sixth  vortex  filaments  (v5  and  v6 
in  Fig.  6)  were  understoodas  near  the  severe  B  VI  locations  and  con- 
vected  downstream  from  the  third  and  fourth  blades,  respectively. 
The  single-trailer  wake  prediction  for  tip  vortex  locations  in  the  top 
view  confirms  the  earlier  2GCHAS  study,'^  '''  and  the  calculation 
consistently  lags  the  measured  data  by  about  two  chord  lengths  for 
an  azimuth  of  35  deg  and  somewhat  less  for  an  azimuth  of  295  deg. 
The  multiple-trailer  wake  model  shows  excellent  prediction  in  the 
top  view,  unlike  the  single-trailer  wake  model. 


The  miss  distance  (side  view)  is  shown  in  Fig.  6  for  the  BL  case. 
The  measured  tip  vortices  cross  the  reference  (passing)  blade.  The 
single-trailer  wake  model  predicts  the  tip  vortex  locations  to  be 
below  the  blade,  hence,  below  the  measurements,  and  the  multiple- 
trailer  model  predicts  the  tip  vortex  to  be  even  lower.  Consequently, 
the  intersections  of  the  blade  and  tip  vortices  are  predicted  to  oc¬ 
cur  more  inboard  than  the  measured  intersections.  Because  of  the 
importance  of  the  miss  distance  to  the  resulting  noise  and  vibra¬ 
tion,  as  well  as  to  the  airloads,  it  will  be  important  to  improve  the 
calculations  of  miss  distance. 

For  the  MN  and  MV  cases,  trends  similar  to  the  BL  case  are 
observed,  as  shown  in  Figs.  7  and  8.  The  multiple-trailerwake  model 
shows  good  correlation  of  the  tip  vortex  geometry  in  the  top  view, 
whereas  the  single-trailer  wake  model  shows  a  lag  of  the  tip  vortex 


Minimum  Noise  (R138) 


Single  Trailer 


X(m) 


Radius  (m) 

a)  Advancing  side,  ^  =  35  deg 


Minimum  Noise  (R138) 
Single  Trailer 


X(m) 


Radius  (m) 


Minimum  Noise  (R138) 


Multiple  Trailers 


X(m) 


Radius  (m) 


Minimum  Noise  (R138) 
Multiple  Trailers 


X(m) 


Radius  (m) 


b)  Retreating  side,  ^  =  295  deg 

Fig.  7  Comparison  of  measnred  and  ealcnlated  tip  vortex  filament  locations  with  single-trailer  and  mnitiple-trailer  wake  models  for  MN  ease. 
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Fig.  8  Comparison  of  the  measnred  and  ealcnlated  tip  vortex  filament  locations  with  single-trailer  and  multiple-trailer  wake  models  for  MV  ease. 


locations  from  the  measured  data.  The  calculated  vertical  position 
is  generally  below  the  measured  position. 

Conclusions 

Calculated  blade  airloads,  motion,  and  wake  geometry  are  com¬ 
pared  with  measurements  from  the  HART-I  rotor.  The  calculations 
were  performed  using  the  comprehensive  analysis  CAMRAD  II, 
with  single-trailerand  multiple-trailer  free- wake  models. 

The  multiple-trailerwake  model  produces  significantly  better  air¬ 
loads  correlation  than  the  single-trailer  wake  model.  It  gives  good 
correlation  for  the  airloads  in  the  BL  case.  In  particular,  the  2-per- 
revolution  normal  force  is  well  predicted,  which  was  not  possible 
in  previous  investigations  using  the  single-trailer  wake  model.  The 


multiple-trailer  wake  model  also  showed  good  correlation  of  the 
measured  and  calculated  wake  geometry  in  the  top  view,  whereas 
the  single-trailer  wake  model  gives  positions  of  the  wake  geometry 
well  forward  of  the  measured  positions. 

However,  not  all  of  the  calculations  were  improved  with  the 
multiple-trailerwake  model.  The  measuredairloads  in  the  minimum 
noise  and  minimum  vibration  cases  exhibit  larger  3 -per- re  volution 
content  than  the  calculated  airloads.  The  measured  and  calculated 
blade  tip  torsion  motion  show  differencesin  2-  and  3 -per- re  volution 
content  that  are  consistent  with  the  differencesin  the  airloads. 

The  miss  distance  predictions  are  deteriorated  with  the  multiple- 
trailer  wake  model  compared  with  the  single-trailer  wake  model 
for  all  BL,  MN,  and  MV  cases.  Even  with  some  poorly  predicted 
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parameters,  it  is  obvious  that  a  multiple-trailer  wake  model  signifi¬ 
cantly  improves  the  understandingof  the  source  of  discrepancy  from 
the  previous  HART-I  correlation  efforts.  Generally  the  calculated  tip 
vortex  positions  are  below  the  measured  positions. 
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